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Abstract 

We motivate the simplest ansatz for the neutrino mass matrix consistent with 
the data from neutrino oscillation experiments, and admitting CP violation. It 
has only two free parameters: an arbitrary mass-scale and a small dimensionless 
ratio. This mass matrix exhibits two symmetries, Democracy and Mutativity, 
which respectively ensure trimaximal mixing of the |z^2) mass eigenstate, and 
mixing parameter values |^23| = 45° and \5\ = 90°, consistent with bimaximal 
mixing of the {u^) mass eigenstate. A third constraint relates the smallness of 
|C/e3p to that of the mass-squared difference ratio, Am^^i/ Am'^^^, yielding the 



prediction sin 6*13 = ^j2Amj^j3Ami^^ ~ 0.13 ± 0.03. 
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1 Democracy and the Solar Data 

The oscihation data from solar neutrino experiments |I| and from the KAMLAND 
experiment [21 all point to a common conclusion, namely that the magnitude, \Ue2\, of 
the relevant MNS lepton mixing matrix element is given by |f/e2p — 1/3- Combined 
with information on atmospheric neutrino oscillations jlj, from the Chooz and Palo 
Verde jH] reactor neutrino experiments and the constraints of unitarity, this strongly 
suggests [Hj that the mixing matrix has one trimaximally mixed column corresponding 
to the 1^2 mass eigenstate: 



^2) = ^ 1 • (1) 




This, in turn, implies that the neutrino mass-squared matrix, := My Ml, in the 
flavour basis^ is an 53 group matrix [Hj (ie. a 3 x 3 'magic square' 0): 
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where a, x, y and z are real parameters and b is complex (we note that is propor- 
tional to the effective Hamiltonian for neutrino oscillations, H = M^/2E, where E is 
the neutrino energy). 

Such a mass matrix, Eq. (j21), is the most general one which commutes with the 
Democracy operator [7]: 

[Ml,V] = (3) 
where the Democracy operator, V, is defined as: 



1 1 1 




(4) 



Hence, the neutrino mass eigenstates (displayed in the Appendix, Eq. ((211)) are all 
eigenstates of V, with |z/2) having eigenvalue V = 1, while {ui) and lz/3) each have 
eigenvalue T> = 0. The Democracy operator generates a group of unitary transforma- 
tions, Uj) = e^"^ [= 1 + P(e*" — 1)] {a a real parameter), under which M^ is invariant: 

UvM^U^ = Ml (5) 

The usual consequences of unitary symmetry follow, eg. the Democracy quantum 
number is conserved in vacuum neutrino oscillations (although, clearly. Democracy 
symmetry is broken by matter effects). 

^Working in the flavour basis, we choose the epsilon phase convention 7 , in which the freedom to 
rephase the charged lepton flavour eigenstates is used to set the imaginary part of the neutrino mass- 
squared matrix proportional to the totally anti-symmetric 3x3 matrix, A, deflned by Aa/s := ea/B-y 
(7^a,/?) a,/3,7= 1,2,3. 
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We note that our approach here is distinct from that in eg. Ref. jH] in which the 
mass matrix itself was taken to be (approximately) proportional to T) and was then 
said to be (approximately) democratic. Here, Democracy refers to the generator of 
the symmetry transformation, and matrices of the (considerably more general) form, 
Eq. (j2)), are precisely invariant under such transformations, Eq. (0). 

The parameters used in Eq. Q are not all observable, in the sense that, for a given 
set of observables, we may always add a given real constant to a and 6, providing 
that we subtract the same constant from x, y and z. This allows us to set Re(6) = 0, 
without loss of generality. Furthermore, only the neutrino mass-squared differences are 
observable in neutrino oscillation experiments, so that the offset of the mass-squared 
eigenvalues is unobservable, and we may set the value of a equal to zero in this context. 
The resulting reduced mass-squared matrix may then be written: 



where d = Im(6) is a measure of CP violation. We note that the two terms 
are respectively pure-real and pure-imaginary. This simple four-parameter ansatz pa- 
rameterises the six independent observables of neutrino oscillations: the solar and 
atmospheric mass-squared differences, Amg^j and Aml^^, the three mixing angles [HI, 
012, 023, 013 and the CP-violating phase, 6. Expressions for each of these are given in 
Eqs. fj3H) - (jnBj) of the Appendix. Two degrees of freedom have been removed by the 
constraint of Democracy invariance, which fixes the magnitudes of the elements of the 
\i>2) mass eigenstate (only two of which are independent), Eq. (^. 

2 Mutativity and the Atmospheric Data 

The data on atmospheric neutrino oscillations j4j, supported by the K2K experiment 
[in] , and from the Chooz and Palo Verde [Sj experiments furthermore indicate that 
the z/3 mass eigenstate is approximately bi-maximally mixed between the /x and r 
flavours. In the three generation context, exact bi-maximal mixing of the z/3 eigenstate 
is understood to mean \U^s\ = \Ut-3\ = l/\/2, or equivalently, |t/^3| = \Ur3\ and 
\Ue3\ = 0. Indeed, taking account of terrestrial matter effects [TT], which tend to 
decouple high-energy electron- neutrino mixing at the atmospheric scale jT2j, it may 
be said that the data on atmospheric neutrinos essentially require jf/^s] — \Ur3\, while 
it is the reactor experiments, Chooz and Palo Verde, which require |f/e3p << 1 (as 
will be addressed in the next section). Taken together, the combined data imply that 
the 1^3 eigenstate is approximately bi-maximally mixed. 

It is by now well-known |1,3] that the first of the two bi-maximal requirements, 
{Ufisl — \Ut3\, can be ensured by a form of /i — r permutation symmetry Tl]. In the 
context of the D-invariant ansatz, Eq. ©, this may be enforced by simply setting 
y = z. It is shown in the Appendix that in the standard parameterisation, this 
corresponds to fixing | sin ^23! = l/"\/2, and |sin5| = 1, exactly what is needed for 



\Ui_i3\ = \Ur3\- We note that for any fixed value of 6*13, this choice corresponds to 
maximal CP-violation (l^asl = 45°, \5\ = 90°). 
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The mass matrix, Eq. © with y = z, now, therefore, has an additional symmetry: 
it commutes with the Mutativity {fi — t reflection) operator: 

[M',-M]=0 (7) 

which is defined as multiplication by the /i — r exchange operator, followed by complex 
conjugation: 




M = CPf,r, where P^^ = 1 (8) 



and C is the complex-conjugation operator. We comment that the refined mass eigen- 
states (those of Eq. (jH)), with y = z, which are given in the Appendix, Eq. (P7j)). as 
well as being eigenstates of Democracy, are now also eigenstates of Mutativity (as 
expected given Eq. with eigenvalues \J^\ = 1. Mutativity being an anti-unitary 
operator, the phases of the eigenvalues are not well-defined, and there is no corre- 
sponding conservation law. There are, however, phenomenological consequences of 
Mutativity symmetry (a symmetry which remains good, even in the presence of matter 
effects). These consequences have been explored thoroughly in a previous publication 
(independent, in fact, of the constraint of Democracy). 

The mixing matrix obtained assuming both Democracy and Mutativity is displayed 
in Eq. ()37|) . and was first introduced in Ref. |15j. It has a single unconstrained mixing 
angle, {Ue^l = sin 6*13 = sin x, in addition to the two observable mass-squared differ- 
ences, AmgQi and Am^^j^. Expressions for these observables in terms of the parameters 
are also given in the Appendix. However, with the assumptions made so far, the z/3 
mass eigenstate is not necessarily bi-maximally mixed (nor even nearly bi-maximally) , 
as jt/esl remains unconstrained, and our mass matrix at this stage represents only an 
intermediate step in our derivation. In fact, this mixing scheme interpolates between 
the earlier tri-bimaximal scheme fHlEI; corresponding to the choice d = in Eq. 
yielding sin6'i3 = 0, and the original trimaximal scheme ^7j, having x = y = z, 
yielding sin6'i3 = l/VS. 

3 Chooz, Palo Verde and the Simplest Neutrino Mass Matrix 

So far, the solar data, primarily, have motivated our introduction of the Democracy 
symmetry, and the atmospheric data, that of the Mutativity symmetry. We still need 
a way to ensure a small value of \Ue3\, as required by the Chooz and Palo Verde 
experiments. We continue in a similar vein, letting the data be our guide. Assuming 
T>- and Al-invariance then, we consider the experimental constraints on the values of 
the three remaining parameters, x, y and d. We find (after inverting Eqs. 
and 
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^ sin' ^ + = — -^atm ( _ ^^Q) 
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d = ^^sin2x, (11) 

where sin^ x = (3/2) sin^ ^^13 < 0.06, as indicated by the results of the Chooz and 
Palo Verde experiments 0, and (Am^^i/Am^^^) = 0.027 ± 0.010 dH. Evidently, 
y (= z) « X, while d < OAx. 

It is clear that the parameter d is at least somewhat smaller than x, and that the 
parameter y (= z) is "very small" in the sense that it is considerably smaller than the 
upper-limit on d. There is no obvious reason why y (= z) should be so small (it follows 
from the smallness of both sin^ x ('^"^soi/'^"^atm)) indicated by experiment). 
We would therefore like to consider the phenomenology of the case in which the very 
small parameter y (= z) is, in fact, exactly zero. This gives six "texture" zeroes jl^ 
in the real part of the (reduced) mass-squared matrix, and results in only two free 
parameters, x and d, remaining: 



= x 1 + -1 1 . (12) 





While the constraint y {= z) = does not obviously increase the symmetry of as a 

whole, it certainly increases that of the real term alone, although it is not clear what 

significance, if any, this has. We will proceed anyway to consider the phenomenology, 

noting that our equations will be at least approximately valid, as long as y{= z) « d. 

By construction, the mass-squared matrix, Eq. ()12j) . still has a single trimaxi- 

mally mixed mass-eigenstate (by Democracy), and \U^3\ = \Ur3\ (by Mutativity), but 

now also builds-in a relationship between the experimentally small mixing parame- 

2 



ter, {Uesl = sin (= y 3 sinx), and the experimentally small mass-squared difference 
ratio, Amgoi/Am^(.j^. In particular, setting y {= z) = in Eq. (fTUj). we see that^: 

sin^ X = ^mlJ Am^^^ (13) 

which yields the new prediction: 



sin ^13 = p/^mlJ?>/\ml^ ~ 0.13 ± 0.03, (14) 

clearly consistent with the constraints from the Chooz and Palo Verde experiments |3] 
(see also eg. [20] for similar, but distinct relations based on different reasoning). The 
smallness of sin ^13 is now ensured and hence, the z/3 eigenstate is nearly bi-maximally 
mixed, as required. 

The observables and parameters are related as follows: 

Am^^i = x±^d^ + a;2 (15) 
Am^^^ = ±2 + x2 (16) 

tan2x = -. (17) 

X 



^This prediction was first made, somewhat obliquely, in 
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For a normal neutrino mass hierarchy, the positive square-root is taken, and x is 
negative. For an inverted hierarchy, take the negative square-root and x positive. The 
MNS lepton mixing matrix is given by Eq. ()37|) . with: 

sin2x = -[l rl, cos^X = -fl + rl- (18) 

From Eq. (|13|). we also have cos^x = Am\2/ ^"fnit^-, where Am32 = Am^^^ — Amg^^j. 
We can now re-write Eq. (|12p succinctly in terms of observables: 

Ml = [cos 2x P,r + ^ sin 2x A] (19) 

where P^r and A are the first and second matrices appearing in Eq. p2|l . respectively. 

We believe this ansatz (or texture) to be the simplest possible one for the neu- 
trino mass-squared matrix, consistent with the data, and admitting CP violation. 
Certainly, the neutrino mass matrix needs an overall scale, and at least one small pa- 
rameter, making a minimum of two parameters (there could be other two-parameter 
schemes, but it is difficult to imagine any with fewer). The six observables of neutrino 
oscillations include 3 large (in fact, ~maximal) mixing angles (6'i2, ^23 and 5), one 
smaller mixing angle (^13 - we know only its upper limit, at present), one very small 
mass-squared difference ratio, Amgoi/Am^t^j, and one overall mass scale, Am^^j^, say. 
These are parameterised here, by a mass matrix, Eq. ()12|) . characterised simply by 
an overall mass-scale, \/d? + x'^ = Am^^j^/2, and a single small ratio, (rf/y^s) ^ 0.2 



(we note that d ^ yAmg^^jAmatm — 0.34a;). The three large mixing angles are thus 
fixed by the Democracy and Mutativity symmetries, while the small and very small 
observables are respectively, ~ linear and ~ quadratic in the small ratio, (^/a/Sx): 

sin^ = {2^mlJ?,^mlJ ^ ]^{d/V?>xf. (20) 

While we cannot say why the ratio {d/ \/?>x) is small, we are left needing to explain 
the smallness of only this single quantity. 

4 Predictions for Long Baseline Neutrino Experiments 

Our expression for sin ^13, Eq. (jlH), enables us to predict appearance and survival prob- 
abilities in future long-baseline neutrino experiments. The leading term in the electron- 
neutrino appearance probability at eg. a neutrino superbeam experiment would be 



P{v,^v:) ^ 2sin2^^i3COs2^i3sin2f^^^^|^^j (21) 

= (0.035 ± 0.013) sin^ (^ ^"^1"^ ^ (22) 

(L is the propagation-distance of the neutrino), while the leading term in the reactor 
anti- neutrino disappearance probability would be: 

1-P{u^^u^) ^ 4sin^gi3Cos^gi3sin^f ^'^^'"^ ) (23) 
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= (0.070 ± 0.026) sin^ [ iE j ^'^^^ 

(it is exactly double -P(z^^ Ve) above, as a consequence of Mutativity symmetry). 
Thus, observable consequeneces are significant, and, in particular, hold out the pos- 
sibility of the measurement of ^^13 in the next generation of long-baseline neutrino 
oscillation experiments j21j . 

Looking further ahead to the prospects for the observation oi CP violation in neu- 
trino oscillations, we note that because the smallness of the solar mass-squared differ- 
ence is linked to the smallness of \Ue2,\ in this scenario, in the limit of CP conservation 
(sin^^i3 — s> 0), Am^Qj also. Hence, we see that the well-known suppression of CP 
violation in neutrino oscillations by two factors, namely the smallness of Amg^j and 
that of l^/esl 122], is attributable to the single small CP- violating parameter [d/\/?)x), 
on which the leading CP-violating term has a ~ cubic dependence. This suppression 
is, of course, somewhat compensated in our present scenario, by the largeness of the 
other three mixing parameters, 612-, O23 and 6, leaving open the possibility to observe 
CP violation in neutrino oscillations at some point in the future. 

5 Discussion 

In this simplest of mixing schemes then, all the features of neutrino oscillation exper- 
iments (aside from the LSND result [22]) are manifest. Both terms of Eq. (|12j) are 
separately symmetric under both the Democracy and Mutativity operators. Taking 
the two terms in order of decreasing magnitude, one may consider the leading term, 
proportional to x, to explain the atmospheric data, its magnitude being directly re- 
lated to the atmospheric scale, Eq. 0, and, on its own, being diagonalised by a simple 
2x2 maximal mixing matrix operating in the /x — r/2 — 3 subspace, thereby ensuring 
sm 

^23 1 = I/V2. Such a mass matrix taken in isolation, however, would result in a 
degeneracy between the z/i and 1^2 mass eigenstates, and zero values for the remaining 
mixing parameters. 

The second term in Eq. ()12j). proportional to d, can be thought of as a perturbation, 
lifting the ui — 1^2 degeneracy with a small eigenvalue difference, Am^Qj ~ d'^/2x, and 
collapsing the i>2 eigenstate into the trimaximally mixed state, Eq. ([Q). It therefore 
explains the solar data, Amg^j, and \Ue2\ = The interplay between the two 

terms in Eq. ()12|) ensures a small sin ^13, of order d/x, allowing CP to be violated. 

We comment finally, that we have not considered here the nature of the neutrino as 
a Majorana or a Dirac particle. Our considerations apply equally well in either case. By 
considering the reduced mass-squared matrix, of Eq. in which we have removed 
the dependence on the parameter a of Eq. (j21), we have ensured no dependence on the 
overall offset of the neutrino masses. The addition to in Eq. (jH)) of any multiple 
of the identity, to represent such an offset (thereby recovering the full M^), clearly 
respects all the symmetries considered here. Our analysis is therefore both insensitive 
to and unpredictive with respect to the phenomenology of neutrinoless double beta 
decay. 
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Appendix: Relationships Between Parameters and Observables 

In order to keep the exposition of the main concepts of this paper concise, we provide 
some of the more complex formulae in this Appendix. 

We start with the most general Democracy-invariant mass-squared matrix, Eq. (j^J. 
In terms of the parameters a, b, x, y and z, the eigenvalues may be written jH]: 



a-Reb- ^3 (Im6)2 + [{x - yf + {y - zf + {z - xf]/2 (25) 

a + 2Reb + x + y + z (26) 

.3 - a - Re 6 + y's (Im6)2 + [(x - y f + {y - zf + (^ - x)2]/2 (27) 

and the MNS lepton mixing matrix is jH]: 



ml 



U 



e 
r 



^1 



3 



V 



V6 



V2 



+ 



V6 ^ V2 

where we have used the abbreviations: Cv 



V2 

1 
1 
1 



i^3 



^x^<P '^^x^4> _i_ ^x^4'~^'^^x^4> 



+ 



V2 

^x^4> '^^x^(t> I 

73 + 



cosx, = smx, C0 



v/6 
COS(/), Sa 



and the phenomenological parameters x ^i-nd are given by: 

VS{z-y) 



tan 20 
tan 2% 



x) + {y — x)' 

\/6 Imb 



[ {x — yY + {y — zY + (-2 — ] 2 



(28) 



sm( 



(29) 
(30) 



The usual mixing observables in the standard (PDG 9J) parameterisation may be 
written in terms of x ^^id 0: 



sin^^ 



13 



sin^^ 



12 



sin^ Q. 



23 



and the CP- violation parameter J 
sin 2x 



J 



cos 2x cos 20) 



1 



3 cos^ 6*13 

1 cos 2x sin 20 

2 ~ 2v^cos2 
is given by: 

Im b 



(31) 
(32) 

(33) 



6^3 3v^[6(Im6)2 + (x-y)2 + (y — zY + (-2 — x) 



(34) 



from which the CP-violating phase 5 can be extracted if required. It is manifest that 
all observables are invariant if we add a given real constant to a and 6, providing that 
we subtract the same constant from x, y and z. Of course all mass-squared differences 
are independent of a. 
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Without loss of generality then, we can choose Re6 = 0, and hence we can write 
the two mass-squared differences: 



^ml,^ -.= 1711 -ml = 2^ d? + \{x - yf + - zf + {z - xf]/2 (35) 

A 2 2 2 



= X + y + z 
+ ^d?+ [{x - yf + - zf + {z- x)2]/2 



(36) 



remain 



where d = a/3 Im6. The expressions for the mixing parameters, Eqs. 
unaltered. 

Invoking Mutativity, ie. setting y = z, Eq. immediately gives tan 20 = (and 
hence sin 20 = 0), so that (we choose sin0 = 0) the mixing matrix is: 



U 



e 
/i 
r 



/ 

V 



1^1 


1^2 


V 3^ 


1 






1 






1 


Ve V2 





2^3 



V2 >/6 
v/2 



(37) 



where 



tan2x 



d 



(38) 



(x -y) 

Substituting in Eqs. we find that this corresponds to standard parameters 

2 



sin^^ 



13 



sin^^ 



12 



sin^ 



23 



We also have: 



sin^ 5 



^"^atm := mg - m 



sm^ X 



3 cos^ 6*13 
1 

2 
1. 



nio — m 



2^d^ + {x-yy 
X + 2y 



+ Jd'^ + ^x-yf 



(39) 
(40) 

(41) 
(42) 

(43) 
(44) 



Finally, setting y{=z) = 0, we arrive at the simplest mass matrix, as discussed in 
the main text, and in particular, Eqs. p5|l - (fT7j) above. 
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